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ABSTRACT 



We study in detail the main features of the unitarized Regge model (CFKS), recently proposed 
to describe the small-Q^ domain. It takes into account a two-component description with two 
types of unitarized contributions: one is the multiple Pomeron exchanges contribution, interact- 
ing with the large dipole size configurations, and the other one consists on a unitarized dipole 
cross section, describing the interaction with the small size dipoles. Its extrapolation to higher 
virtualities is performed, analyzing the ratio between soft and hard pieces and comparing the 
resulting dipole cross section to that from the saturation model. Diffraction dissociation is also 
considered, showing the scaling violations in diffractive DIS and estimating the corresponding 
logarithmic slope. 
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1 Introduction 



The study of a new regime of QCD, that of high density of partons, has drawn much attention 
in the last years. The key discovery was the observation at HERA of the fast growth of parton 
densities (mainly gluons) as the energy increases in experiments of deep inelastic scattering. 
Taking (t*°* ~ (F2 ~ a;""(°)+^), values of A = a(0) - 1 in the range 0.1 - 0.5 have been 

reported, depending on the virtuahty of the photon. However, this steep growth should be 
tamed, leading to the expected limit given by the Froissart bound (a < (logs)^ as s — >• cxD)[jl|. 
This boundary has been derived from very general properties of the S-matrix, namely unitarity. 
A cross section growing as any positive power of s would violate unitarity at asymptotic energies. 
Thus, theoretically, some kind of saturation of this growing due to unitarity effects is expected 
. The dynamics of such very dense partonic systems is very interesting and has been studied 
by many authors both in DIS Q and in high energy nuclear interactions 

The description of the 7*p collision in the frame where the proton is at rest is very appro- 
priate to include unitarity corrections. In this frame, the virtual photon 7* emitted by the 
incoming lepton fluctuates into a qq pair. This system then suffers multiple interactions with 
the proton. Such multiple interactions restore unitarity even in the case where it would be 
violated in a single collision. In the model developed in [0, ^, all these corrections have been 
taken into account, and their strength is constrained by diffractive data. Therefore, the ratio 
(jdtff jiytot jg rgia^gfj unitarity corrections. This is a common feature to any realization of 
the Gribov model [0, where the amount of rescatterings is related to diffractive production by 
means of AGK-cutting rules 0. 

In parton language, the increasing number of gluons in a proton as x — >■ makes gluon 
fusion very probable. This fusion produces gluons of higher longitudinal momentum, stopping 
the growing of those with the smallest x. In this way unitarity is not violated. Such a procedure 
was implemented on theoretical grounds from QCD through the multiladder exchange using the 
GLR formalism 0, giving rise to non-linear effects in the standard linear DGLAP approach. 
The outstanding quantity emerging from the unitarization procedure is the saturation scale 
X, 6), setting the region where saturation phenomenon starts to be meaningful. The 
QCD-inspired phenomenological model 0, for instance, introduces a quite clear identification 
for this scale Ql{x) ~ 1/Rq{x). There the saturation radius Rq{x), related with the mean 
transverse distance between partons, is properly extracted from the small-x data from HERA. 

In any of the descriptions, the unitarity corrections are given by the non-linear terms, 
and a phenomenon of saturation is expected when these terms become important. Since the 
gluons are the partons driving of the high energy processes, the signals of the saturation effects 
should appear in the observables probing the gluonic content of the proton (or the nucleus) 
[Pm]. In the nuclear case, the gluon density is ~ A^^^ higher than in the proton. This makes 



unitarity corrections more important for nuclei, producing the well-known shadowing of F2 |TT 
Saturation will thus start at smaller energies in nuclei than in protons. Such a fact is the main 
reason for the increasing interest in the forthcoming eA experiments, where the nucleus will be 



studied at energies higher than currently available 
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The open question is if the unitarity corrections have aheady shown up at present energies 
and if the saturation has been reached. In particular, at HERA, they should appear in the 
small-a; and small-Q^ data |jl2| . There are several proposals in this direction , mainly for 



the case of heavy-ion collisions |]T3, but a definitive answer is still missing. The main difficulty 



that we are faced with is the saturation scale Qs, staying in the transition interval of 1-2 
GeV, which leads the effects to be hidden in more inclusive observables. In this kinematical 
region the standard QCD perturbative expansion is expected not to be completely reliable. 
For instance, higher twist terms to the linear approach should be taken into account in such 
a domain. Moreover, this is the transition region between the soft and hard domains, i.e. 
the perturbative approaches (including saturation or properly adjusting initial conditions) and 
the Regge- inspired models are competing, and both frameworks seem to describe the current 
small- a; data. 

Bearing in mind that the saturation phenomenon is required in a complete understanding 
of the high energy reactions, and that a consistent treatment of both inclusive and diffractive 
processes should be taken into account, in this work we study derivative quantities using the 
Regge unitarized CFKS model ^. In this hybrid model, both soft (multiperipheral Pomeron 
and reggeon exchanges) and hard (dipole picture) contributions are properly unitarized in an 
eikonal way. This approach describes the transition region and can be used as initial condition 
for a QCD evolution at high virtualities The extrapolation to the higher-Q^ domain is 



also performed here, checking the behaviour of the model without including QCD evolution. 
We discuss the similarities and/or connections with the phenomenological saturation model 
0, stressing that a QCD evolution is required for a correct description of higher in the 
inclusive case. For the diffractive case, such a procedure is not formally required, since the non- 
perturbative sector is dominant in this case. The diffractive structure function is extrapolated 
to the available larger-Q^ range. In particular, the diffractive logarithmic slope, which has 
been claimed as a possible new observable to disentangle dynamics ||16|, [l^, is calculated and 
compared with the result from the saturation model. 



2 The inclusive case 



We start by briefly reviewing the CFKS approach. It interpolates between low and high virtual- 
ities Q^, which are related to the dipole separation size, r, at the target rest frame, considering 
a two-component model P, |^ . Considering the unifying picture of the color dipoles, the sepa- 
ration into a large size (in |^ it is called L) and a small size (called S* in [Q) components of the 
qq pair is made in terms of the transverse distance r between q and q. The border value, ro, 
is treated as a free parameter - which turn out to be tq ~ 0.2 fm. This value agrees with the 
correlation length of nonperturbative interactions observed in lattice calculations. Hereafter we 
use the notation soft for the large size conflguration and hard for the small size one. This sep- 
aration corresponds, in Ref . , to the separation of the qq pair fluctuation into the "aligned" 
component (large transverse size, soft component) with a strongly asymmetric sharing of the 
momentum fraction between q and g, and the "symmetric" component (small transverse size, 
hard component). 
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The soft component considers multiple Pomeron exchanges (and reggeon /) implemented in 
a quasi-eikonal approach |]18|. It also includes the resummation of triple Pomeron branchings 



(the so-called fan diagrams). The initial input is a phenomenological Pomeron with fixed 
intercept ap(0) = 1 + ejp = 1.2 (further changes are due to absorptive corrections), and 
an exponential parametrization for the t dependence is considered. In the impact parameter 
representation, the 6-space, it looks like (in photoproduction Q"^ = 0): 

X'^{s,b) ^ CjpJ^ f-V" eM-b'/B^^is)] , (1) 



Bel{S) \So 

where Bgi{s) is the elastic slope, which is parametrized as in the hadronic reactions. The /jp 
is an effective Pomeron-proton coupling. In the electroproduction case, the initial input is 
described in an analogous way: 

corresponding to the Regge parametrization for the amplitude of the soft Pomeron exchange. 



similar to the Donnachie-Landshoff one [|T9|. The function i?(x, Q^) comes from the exponen- 
tial assumption about the t dependence and further transformation to the impact parameter 
representation. We shall remark here that, in the CFKS approach, the authors consider a 
Pomeron fixed intercept of ap(0) = 1 + ep = 1.2, ejp = 0.2 (a semi-hard value rather than a 
soft one). The unitarization effects, described by multi-Pomeron exchanges, lead to an effective 
intercept e^ff = '^^"■f^n? ^ that decreases as or x decreases due to the increase of shadowing 
effects. 

The resummation of the triple-Pomeron branches is encoded in the denominator of the 
amplitude the Born term in the eikonal expansion. Moreover, the corrected amplitude 

is eikonalized in the total cross section, 

x"^(x,Q^6) = x'^{x,Q'^,b) 

1 + ax3(a;,Q^&) ' 



a^^{x,Q',b) l-exp[-x"^(x,Q^6)J . (4) 

where the constant a depends on the proton- Pomeron and the triple-Pomeron couplings at zero 
momentum transfer {t = 0). Refs. P, give a more detailed discussion. 

The eikonalization procedure modifies the growth of the total cross section from a steep 
power-like behavior to a milder logarithmic increase. The above parametrization corresponds 
to the interaction with the large size dipole configurations and therefore dominates in low-Q^ 
values. The total soft contribution is obtained by integrating over the impact parameter the 
cross section at fixed b, a"'^{x,Q'^,b), 

a'°f\s, Q')=aI d% a'°f\s, Q\ b) . (5) 



The hard component is considered in the color dipole picture of DIS |2D|. The dipole cross 
section, modeling the interaction between the qq pair and the proton, o"'^*^°'^(x, r), is taken from 
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the eikonalization of the expression above x"^(s,6, Q^) aheady corrected by triple-Pomeron 
branching (the fan diagrams contributions). The configurations considered are those with a 
small transverse distance between the quark-ant iquark pair in the dipole. The corresponding 
cross section is extracted by considering the contributions coming from distances between 
and ro = 0.2 fm (1 GeV^^) - see discussion above whereas for r > ro the contributions are 
described by the soft piece already discussed. In such small distances, perturbative QCD is 
expected to work. The total cross section considering this dipole cross section is expressed as 
i: 

<yt:''M') = I <fr j^da\^^/^{a,r)\'a'l^f^%{x,r), (6) 
^Itsi^.r) = 4 Jd%a-^{x,Q',b,r), (7) 
a"^(x,g^6,r) ^ l-exp[-rV^(a:,g^&)], (8) 

where T and L correspond to transverse and longitudinal polarizations of a virtual photon, 
\E^i^iq(a,r) are the corresponding wave functions of the gg-pair. 

The dependence is introduced in the Born term of the eikonal expansion, presented in 
the last expression above, in order to ensure the correct behavior determined by the color 
transparency: for small r the growth in radius should be proportional to r^, a"'^{x, Q^, b, r) ~ 
r'^ f{x,Q'^,b). This condition, valid for fixed s and as r ^ 0, is a property of the single 
Pomeron exchange. Thus a factor has been introduced in the eikonal of eq. (8). 

Another difference between the soft and hard component is the fact that the contribution 
of the /-exchange (reggeon exchanges) to the hard component is very small, and has been 
neglected. 

The weight of each contribution (soft and hard) in the total cross section [and F2{x,Q^)] 
can be obtained, providing an analysis of the role played by each piece of the model. Such a 
procedure allows us to explicit the regions of x and where the sectors contribute. In Figs. 1 
and 2 we calculate the ratio Rsoft, defining the fraction of the total contribution arising from 
the soft sector: 



RsofAx, g2) = ^"'^ ^ . (9) 

From Fig. 1 we note that the soft contribution slowly increases as the momentum fraction 
X goes to higher values, almost independently of the virtuality Q^. This is due to the fact that 
higher reggeon trajectories / are included in the soft part, but not in the hard one. Regarding 
fixed virtualities, the soft piece dominates completely the total cross section at = 0.045. As 
Q'^ increases the contribution goes down. For instance, at = 10 GeV^ it contributes about 
half of the cross section. Extrapolating up to higher virtualities, the soft piece saturates at 
about 5-15% of the total result. 
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Figure 1: The ratio Rsoft as a function of x at fixed virtualities. 



Figure 2 clearly shows that the soft piece is dominant at = 0.01 and decreases as the 
virtuality grows. The behavior is monotonic, almost independent of the momentum fraction x. 
For instance, at Q"^ = 100 GeV^, it contributes with 20% at x = 10~^ and 5% at x = 10~^. 
Such a reduction on the soft content is related to the coupling of the photon to the asymmetric 
dipoles gloft^Q"^) ~ V(l + Q'^ I ''^loft) and to the enhancement in Q"^ provided by the photon 
wave function (at high ^ Qli^) the symmetric dipole configuration provides the scaling 
with logarithmic violation). 

An interesting issue is the relation between the dipole cross section coming from the CFKS 
model and the phenomenological one of G.-Biernat-Wiisthoff 0. The GBW cross section is 
parametrized as: 



a 



GBW t 



X, r) 



Rlix) 



X 



1 - exp(-rV4i?o(a;)) 
-VgcV-^ 



(10) 
(11) 



where ctq = 23.03 mb properly normalizes the dipole cross section. The remaining parameters 
are A = 0.288 and xq = 3.04 x 10~^, all of them determined from the small- a; HERA data. The 
Ro{x) is the main theoretical contribution, defining the saturation scale, which is related with 
the taming of the gluon distribution at small x (unitarity effects) . The above expression has 
been used to describe both inclusive and diffractive structure functions, in good agreement with 
the experimental results. The comparison between this approach and the CFKS dipole cross 
section is shown in Fig. 3. We have plotted the adimensional result, since the normalization 
for the CFKS dipole cross section, o"o, is not determined from data. Indeed, for a comparison 
with experiment using only the hard piece from CFKS, the adjustable parameters would have 
to be refitted. We consider here that this can be absorbed by a suitable normalization, and 
carry the r interval beyond the range set by the model (r < tq). The main feature of the GBW 



5 



x=10 



L 
10 ' 



10" 



10" 



10' 



10^ 



Figure 2: The ratio Rsoft as a function of at fixed momentum fraction x. 

parametrization is that it ensures that the dipole cross section grows hnearly with at small 
transverse separation, whereas it saturates at large size configurations. The picture emerging 
from the CFKS is slightly different, presenting a mild (logarithmic) increase with r, away from 
huge separation sizes that shifts the saturation scale up to very high virtualities. Although 
the continuous and smooth increasing with the radius, in the CFKS approach the cross section 
underestimates the GBW one for all r. 

A comment on the normalization is in order. The GBW formula would correspond to 
the hard part of CFKS without triple-pomeron (a = 0), and taking a step function for the 
profile instead of a gaussian. This makes unitarity corrections stronger. In any case, we 
can compare GBW and the hard part of CFKS: taking a = and exp[— Q^)] — 
0[6^ — R^x.Q"^)] in CFKS, doing the integral in h and comparing with GBW, one obtains 
(Jcf'ks{^^ r) = (To 1 — exp(— X^{s^ b, Q^) , with ctq = uRi^x, Q^) ~ 20 mb, in agreement with 
GBW value. This value, however, depends logarithmically on x and Q^, because of the increase 
of the proton radius, which is taken into account in (2). The comparison between the exponents 
(the eikonals, which contain most of the parameters in CFKS case) of the above expression and 
of the eq. (10) is less clear, as the x— dependences are different and the triple pomeron can not 
be neglected in this case. 



3 The diffractive case 

The diffractive sector in the CFKS approach is constructed by a three-component model f^, ^ , 
using the AGK cutting rules to relate the elastic multiple scattering amplitude to the inelastic 
diffractive contribution ||^. The first term comes directly from the soft piece, the second one 
from the triple-Pomeron (and the reggeon /) interaction and the last one from the hard (dipole) 
piece. We notice that these contributions define only the energy, s (and momentum fraction 
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Figure 3: The comparison between the saturation dipole cross section from G.-Biernat- 
Wusthoff (GBW) and CFKS as a function of the transverse dipole separation r at fixed x 
is). 

x), and the virtuahty dependences. The spectrum on (3 is introduced by hand, based on earher 
soft and hard (pQCD) calculations. The first component is written as: 



F.lomix, Q^ f3) - F^Jf-^ K,{s, Q') (1 - /?)"^(«^) , (12) 

where (Born)^^^^ ~ X"'^{s,Q'^)x"'^{s,Q^) is the lowest-order (Born) approximation for that 
function, with i,k = IP, f. The suppression factor due to higher order multipomeron exchanges 
is K^is, Q2) = a;°}f with a% = AgHQ^) J d% [a,oft{s, Q\ h)f. Further details can be 
found in 0. The (3 dependence is taken from the typical CKMT Pomeron structure function. 



which is connected with the deuteron structure function by the identification x — ^ /5 
The hard contribution is expressed as: 

F?iHar,){^. Q\ /?) ~ ^yfar, /^'(l " 2/5)^ + ofj, i3\\ - (5) , (13) 



where the (3 dependence comes from a pQCD guess for the Pomeron structure function p2 |. 
Also, (jfl^d^ = J d% [o"^a^rf(s, Q^, hyf'. However, the (3 spectrum is slightly different from the 
most recent pQCD calculations, where the transverse contribution behaves like ~ 13{1 — j3) 



2^]. The hard component dominates at large (3 in the CFKS approach, mainly concerning 
the transverse component, where the charge in the j3 dependence will no significantly modify 
the data description, moreover due to few data to constraint the adjustable parameters in this 
region. 

Regarding the /3 dependence, the region for medium values (/3 ~ 0.4) is dominated by the 
soft term, which in pQCD is associated to the transverse photon contribution The small 
j3 region is dominated by the triple-Pomeron piece, in agreement with the pQCD expectations. 
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Figure 4: The diffractive structure function as a function of at fixed Mx and W . The 
prehminary data are from ZEUS Collaboration (triangles). The published data are the circles. 
The CFKS and GBW results are shown in the same plot. 

which is obtained by considering the higher twist gg+gluon configuration. Moreover, the hard 
contribution is leading in the large-/? region, associated in this case with a suppression of the 
transverse contribution and an enhancement of the longitudinal piece in comparison with the 
expected pQCD behavior |23]. 



The CFKS approach describes with good agreement the diffractive DIS data in the broad 
range < < 18 GeV^. In order to study the model in comparison with the pQCD ap- 
proaches, here we extrapolate the prediction for the diffractive structure function of the CFKS 
approach for higher values of the virtuality. We use the preliminary ZEUS analyses, considering 
the dependence at fixed mass Mx and centre-of-mass energy W These data provide 
information at both small and large virtualities bins. It is interesting to compare the predictions 
of the CFKS model and the saturation model P] for the diffractive structure function. Both 
models are depicted in the plots of Fig. ^. 

The agreement of CFKS approach with data is remarkable even at high virtualities, where 
the model is expected not to be reliable. However the interpretations at low are quite 
different. In the saturation model, the reliability of the pQCD calculation is extended to 
smaller virtualities through the saturation scale Rq{x)^ where the dependence is mostly due to 
the longitudinal photon configuration, by the higher twist gg-|-gluon. Instead, in the CFKS 
model the main contribution in the region of interest comes from the soft triple-Pomeron 
contribution. 

As a final study, we perform the calculation of the logarithmic slope of the diffractive 
structure function F^^^\ The motivation is that this observable is a potential quantity to dis- 
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tinguish soft and hard dynamics in diffractive DIS |16, 17] . The measurements of the derivative 
quantity F2-slope on have allowed a renewed interest on testing the matching of hard and 
soft approaches and have provided constraints for the saturation formalisms. The reported 
turnover on the x dependence has been associated with the transition region between the soft 
and hard domains. When we focus in diffractive DIS, in particular the structure function F2 ■, 
the situation is far from clear: initially considered as a predominantly soft process, the ex- 
perimental results suggest that the diffractive cross section at HERA contains hard and soft 
components. 

Still, the diffraction stands a more profitable field to study saturation effects than the 
inclusive case. This comes from the fact that, in DDIS, the large dipole size configuration 
(soft content) is more relevant than in the DIS reaction^. Although there are quite different 
approaches, based on different physical dynamics, applied to the interpretation of the diffractive 
measurements, almost all of them fit the data set properly ||25||. Therefore a derivative quantity, 
the diffractive logarithmic slope, has been proposed. It can help to distinguish the underlying 
dynamics in diffractive DIS, settling the validity range of the different approaches, if such 
observable is measured. Here, we have calculated the slope as a function of the Pomeron 
momentum fraction and we have performed a comparison between the CFKS and the G.- 
Biernat-Wusthoff (GBW)approaches. 

In this calculation, the main feature of the GBW approach is the presence of positive and 
negative slopes, contrary to the pQCD non-saturated case |T^. In the pQCD model without 
saturation, the parameters fixed by the previously available data lead to a predominantly 
positive slope for all kinematical region, even at large (large We notice that this situation 



can be changed by a further analysis, considering the updated measurements of HI p5|, which 
have enlarged the available kinematical range and have provided new measurements for the 
regions previously covered. 

It is important to emphasize that a transition in the slope takes place in the preliminary 
ZEUS analyses of diffractive DIS at large virtualities, where the saturation model [§] is 
considered to describe the dependence of the diffractive structure function. In this model, 
the analysis is performed using Mx and W as kinematical variables instead of /3 and xjp, due to 
the similarity of the behaviors of dajdMx and o"toi(7*p) in the same kinematical range. There, 
the growth of xipF-f versus is stopped at ~ 10 GeV^ and decreases smoothly for larger 
virtualities. The transition region corresponds to /3 ~ 0.2 for Mx = 5 GeV and (3 ~ 0.07 for 
Mx = 11 GeV [l^], /3 = QV(Q^ + ^D- These features can be verified observing the plots of 
Fig. (4) at large virtualities. 

Here, we compared the results from the saturation model with those from the model analyzed 
in this work. They are shown in Fig. (|^), where the slope is calculated as a function of xjp for 
fixed (3 and at two different virtualities. We choose = 1 and 10 GeV^ because this is the 
region where the CFKS model is formally valid. However, we emphasize that it can be extended 
to higher virtualities in the diffractive case, since the soft component is stronger than in the 
inclusive case. The saturation model produces a transition between positive and negative slope 
values at low P = 0.04 (upper plots), while it shows a positive slope for medium and large /?. 
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Figure 5: The logarithmic slope of the diffractive structure function as a function of a;jp at 
fixed f5, for Q'^ = 1 and Q"^ = 10 GeV^. The CFKS and GBW results are shown in the same 
plot. 



Instead, the CFKS approach presents a positive slope for the whole and xjp ranges, flattening 
at large /3, similarly as the non-saturated pQCD calculations. These features can be probed if 
the slope is measured and could help the understanding of the underlying dynamics. Finally, 
the results above, mainly those ones for the structure function, corroborate the CFKS 
model as a consistent hybrid approach to describe diffractive DIS, with a close connection to 
the above mentioned pQCD calculations. 



4 Conclusions 

A deeper understanding of the saturation phenomenon is required to perform reliable estima- 
tions for the current and forthcoming high energy reactions. The saturation scale, which sets 
the onset of the unitarity corrections, is found to be in the transition regime of low x and Q^. 
In this domain, both Regge-inspired phenomenology and improved pQCD calculations (pertur- 
bative shadowing, higher twist), considering unitarity effects, are able to describe the data with 
high precision. The most advantageous ones are those describing in an unifled way the inclu- 
sive processes as well as diffractive ones. In this letter we have considered the two-component 
multireggeon model of p, and calculated some related quantities. 

The ratio of the soft content in the model has been calculated, verifying that it dominates 
at low Q^, diminishing at higher virtualities. This shows that the unitarity corrections in 
this model are more important in the soft component than in the hard one. Moreover, these 
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corrections are higher twist at large Q"^ in the second case. 



We have also studied the robustness of the CFKS model to describe a large range in 
without to consider a pQCD evolution. A good hint to answer this question comes from the 
analysis of the hard piece (symmetric dipole configurations), in particular the corresponding 
dipole cross section. We have found that the r-saturation of this quantity lies larger at values of 
the radius than in the phenomenological GBW model. A similar consideration is far from clear 
for the diffractive case, where the non-perturbative (soft) component plays a more important 
role. 

We have extrapolated the estimations for the diffractive structure function at high virtuali- 
ties, verifying that a broad description is obtained, and that it is in reasonable agreement with 
the saturation pQCD model. An additional quantity has been proposed in order to describe 



the dynamics of the diffractive dissociation |ll6| , p!7|] , in particular the diffractive slope. It has 



been calculated using the CFKS model, and its main feature is a similar behavior to the one 



predicted by pQCD calculations |23 . 
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